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C 3 1 A strongly momentum-dependent dressed-quark mass function is basic to QCD. It is cen- 

tral to the appearance of a constituent-quark mass-scale and an existential prerequisite 
for Goldstone modes. Dyson-Schwinger equation (DSEs) studies have long emphasised 
this importance, and have proved that QCD's Goldstone modes are the only pseudoscalar 
C 3 ■ mesons to possess a nonzero leptonic decay constant in the chiral limit when chiral sym- 

metry is dynamically broken, while the decay constants of their radial excitations vanish. 
Such features are readily illustrated using a rainbow-ladder truncation of the DSEs. In 
this connection we find (in GeV): f Vc (is) = 0.233, »7l^ c (2S) = 3.42; and support for in- 
terpreting ?7(1295), r;(1470) as the first radial excitations of 77(548), 77' (958), respectively, 
' and X(1460) as the first radial excitation of the kaon. Moreover, such radial excitations 

have electromagnetic diameters greater than 2 fm. This exceeds the spatial length of 
lattices used typically in contemporary lattice-QCD. 
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The meson spectrum contains three pseudoscalars [I G (J P )L = 1~(0 _ )5] with 
masses below 2 GeV: 7r(140); 7r(1300); and 7r(1800). In the context of a model 
constituent-quark Hamiltonian, these mesons are often viewed as the first three 
members of a QQ n^So trajectory, where n is the principal quantum number; i.e., 
the 7r(140) is the S'-wave ground state while the others are its first two radial ex- 
citations. This reasoning suggests that the properties of the 7r(1300) and 7r(1800) 
are likely to be sensitive to details of the long range part of the quark-quark in- 
teraction because the constituent-quark wave functions will possess material sup- 
port at large interquark separation. Hence the development of an understanding of 
their properties may provide information about light-quark confinement. The pseu- 
doscalar trajectory is particularly interesting because its lowest mass member is 
QCD's Goldstone mode. Therefore an explanation should simultaneously describe: 
(1) chiral symmetry and its dynamical breaking; and (2) the possible correlation 
of the trajectory's higher mass members via an approximately linear radial Regge 
trajectory. Outcome (2) does not require that confinement in light-quark systems 
be expressed through the formation of a flux tubePlt is easily obtained in Poincare 
invariant quantum mechanic^ but requirement (1) is not. 
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A Poincare covariant and symmetry preserving treatment of quark-antiquark 
bound states can be based on the Bethe-Salpeter equation (BSE) 

,A 

T tu (k;P)= [ X (q;P)] sr K%(q,k;P), (1) 

Jq 

where: k is the relative and P the total momentum of the constituents; r,. . . , it repre- 
sent colour, Dirac and flavour indices; x(q; P) :— S(q + )T(q; P)S(q-), q± = q±P/2; 
and J represents a translationally invariant regularisation of the four-dimensional 
integral, with A the regularisation mass-scale. In Eq. JU, S is the renormalised 
dressed-quark propagator and K is the fully amputated dressed-quark-antiquark 
scattering kernel. 

The dressed-quark propagator appearing in the BSE's kernel is determined by 
the renormalised gap equation (the quark Dyson-Schwinger equation [DSE]) 

S(p)- 1 = Z 2 (» 7 • p + m bm ) + Z 1 J\ 2 D^(P ~ ?)y7,%)r;(9,p) , (2) 

where is the dressed-gluon propagator, Y v (q,p) is the dressed-quark-gluon ver- 
tex, m bm is the A-dependent current-quark bare mass, and Z\p{(, 2 , A 2 ) are, respec- 
tively, the quark-gluon- vertex and quark wave function renormalisation constants. 
The gap equation's solution has the form 

Sip)- 1 = [z7-p + M(p 2 )]/Z(p 2 ,C 2 )- (3) 
It is obtained from Eq. J5J) augmented by the condition S(p)~ 1 \ p2 _ ( * 2 = ■ p + 
m((), where m(Q is the renormalised mass: Z 2 (( 2 , A 2 ) m bm (A) = Zi(( 2 , A 2 ) m(£) , 
with Z4 the Lagrangian mass renormalisation constant. The chiral limit in QCD is 
unambiguously defined by51z 2 (C 2 ,A 2 )m bm (A) = 0,VA 3> C: which is equivalent to 
stating that the renormalisation-point-invariant current-quark mass rh = 0. 

It is a longstanding prediction of DSE studies in QCD that the dressed- 
quark propagator receives strong momentum-dependent corrections at infrared 
momenta! ^ ' v l vl This phenomenon is illustrated in Fig.[U The prediction was con- 
firmed in recent simulations of lattice regularised QCE and the conditions under 
which pointwise agreement may be obtained have been explorediSECQlIII 

The pion's properties are fundamentally governed by the phenomenon of dy- 
namical chiral symmetry breaking (DCSB). One expression of the chiral properties 
of QCD is the axial-vector Ward-Takahashi identity 

P,T^(k;P) = S-\k+)i l5 ^ + ?75 ^5- 1 (fc_) - 2im(0rUk;P), (4) 

which is here written for two quark flavours, each with the same current-quark mass: 
{t 1 : i — 1,2,3} are flavour Pauli matrices. In Eq. @, T 3 5 (k; P) and T j 5 (k; P) are, 
respectively, the renormalised dressed-axial-vector and -pseudoscalar vertices, each 
of which satisfies a unique inhomogeneous extension of Eq. £[} . Equation Q is an 
exact statement about chiral symmetry and the pattern in which it is broken. Hence 
it must always be satisfied. Since that cannot be achieved veraciously through fine 



February 8, 2008 23:58 WSPC/INSTRUCTION FILE cdroberts 



On the complexion of pseudoscalar mesons 3 



10 l 



% 10 



Q. . 2 



10" 



10 



















— b 






■— c 
-- s 


\ \ 

S X. 




— u,d 

- - chiral limit 


\ ~z 
\ ' _ Z 






\ 






\ 

s 






N 

s 















10"' 



10" 



10 u 
p 2 (GeV 2 ) 



10' 



10' 



Fig. 1. Quark mass function obtained with the quenched gap equation kernel described in Ref. 
At a renormalisation point £r = 19 GeV, the current-quark mass values are: m u (j(Ck) = 3.7 MeV, 
m s (Cfl) = 82 MeV, m c (( R ) = 0.97 GeV, and m b (( R ) = 4.1 GeV. The indicated solutions of p 2 = 
M(p 2 ) 2 define the Euclidean constituent-quark mass (in GeV): M^=o = 0.27, M U)< j = 0.28, 
M s = 0.36, M c = 1.3, and Mf, = 4.5. In perturbation theory M^— o(p 2 ) ^ is impossible in the 
chiral limit. Its appearance here is innate to dynamical chiral symmetry breaking. The ultraviolet 
behaviour of the chiral-limit mass function is qualitatively different from that of the solutions 
obtained with a nonzero current-quark mass For p 2 < 1 GeV 2 the u, d mass function is almost 
indistinguishable from the chiral-limit result: this is the essentially nonperturbative domain. 



tuning, the distinct kernels of the DSEs involved; namely, Eq. and the relevant 
inhomogeneous extensions of Eq. , must be intimately related. Any theoretical 
tool employed in calculating properties of the pseudoscalar and pseudovector chan- 
nels must preserve that relationship if the results are to bear fidelity to QCD. A 
weak coupling expansion of the DSEs yields perturbation theory and satisfies this 
constraint. However, that truncation scheme is not useful in the study of intrinsi- 
cally nonperturbative phenomena. Fortunately a nonperturbative systematic and 
symmetry preserving scheme exists. (Reference gives details.) This entails that 
the full implications of Eq. Q can be both elucidated and illustrated. 

Since every flavour nonsinglet pseudoscalar meson is exhibited as a pole contri- 
bution to the axial- vector and pseudoscalar vertices, it follows from Eq. 0} that the 
residues at these poles satisfy-"^ 

A B n4 B = 2m(0/V B (C). (5) 

wherein "n" is the principal quantum number of the n^o trajectory with, e.g., 
n = corresponding to the ground state 7r(140) pseudoscalar meson, and 

U n & = Z 2 tr | A V 75 7„ Xl„ fa P ), ip*J£) 6 ij = Z, tr J i r S xl n («5 P) , 

(6) 
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both of which are gauge invariant and cutoff independe nt. T he identity expressed 
in Eq. is exact in QCD, and valid: for every 0~ mesonpland irrespective of the 
magnitude of the current-quark mass.^l 

For the ground state meson, DCSB entails /° := lim^^o fir and p° (C) := 
limm^o p(0 — ~(QQ)(/ f^o 0- Hence the Gell-Mann-Oakes-Renner relation is 
clearly a corollary of Eq. JSJ) and the ground state pion is QCD's Goldstone mode.EH 

For the n > 1 pseudoscalar mesons: rn Vn>1 > m no by assumption, and hence 
m^ n>a 7^ in the chiral limit; the existence of a bound state entails that Xn(k; P) is 
a finite matrix-valued function; and, furthermore, the ultraviolet behaviour of the 
quark-antiquark scattering kernel in QCD guarantees that 

pl n {0-= lim n ^„(C)<oo, Vn. (7) 

m — >(1 

Hence, it is a necessary co nseq uence of chiral symmetry and its dynamical breaking 
in QCD; viz., Eq. ©, thaJSI 

/2 B = 0,Vn>l. (8) 

These outcomes are legitimate in any theory that has a valid chiral limit. It 
is logically possible that such a theory does not exhibit DCSB; i.e., realises chiral 
symmetry in the Wigner mode. This is the circum stan ce, e.g., in QCD above the 
critical temperature for chiral symmetry restoration!^ Equation JSJ) is still valid in 
the Wigner phase but its implications are different!^ In this event cx m; i.e., 
the leptonic decay constant of the ground state pseudoscalar meson also vanishes 
in the chiral limit. Moreover, since f wo Ptt (C) "<x — (<7<z)° an d (qq)® W oc m 3 in the 
Wigner phasepl it follows that if a rigorously chirally symmetric theory possesses 
a massless pseudoscalar bound state then m^ Q "cx rh. In this case there is also a 
degenerate scalar meson partner whose mass behaves in the same manner. 

The exact results just reviewed can be illustrated accurately in a model that 
both preserves QCD's ultraviolet properties and exhibits DCSB. One such is the 
renormalisation-group- improved rainbow-ladder DSE model presented in Ref. |17j . 
whose widespread application is reviewed in Ref. ^B] ■ The utilisation of this model 
as an illustrative tool in the present context is described in Ref. ^1] and Table ^ 
epitomises the calculations. In regarding the table one should note that the rainbow- 
ladder truncation gives "ideal" flavour mixing and the study used rh u = fhd ^ m s ^ 
rh c . Thus the second row simultaneously describes four degenerate pseudoscalar 
mesons for each n, namely, a {ud, ud — du,du\ isotriplet and a ud + du isosinglet; 
the third row describes ss mesons; and the fourth, cc mesons. The isotriplet mesons 
of Row 2 may plainly be associated with the n 1 ^ pion- like radial trajectory but 
what should one think of the isosinglet, and what of the ss composites of Row 3? 

Th ese q uestions are answered by a nonzero value of the topological susceptibility 
m Qcrjnibut this feature is not realised in the rainbow-ladder truncation. However, 
starting from the model of Ref. [22] one mav easily explore the resolution of the U(l) 
problem via a Witten-Veneziano-likc mechanism; viz., one describes a rudimentary 
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Table 1. Results for the - ground and first radially excited states: m(£o) := "i/(ln£o/AQCD) 7m , 
7 m = 12/25, Co = 1 GeV. As a guide, we list experimental values reported in Ref. 1191 (in 
GeV): m^ = 0.14, m ni = 1.3 ±0.1; m v = 0.548, = 0.958, m^ (2S) = 1.294 ± 0.004; 

m v'(2S) = 1.476 ± 0.004; m Vc(ls) = 2.9796 ± 0.0012, m Vc(2S) = 3.65 ± 0.01; /„■„ = 0.092; andUl 
f Vc = 0.237 ± 0.053. Dimensions are GeV except for p((), GeV 2 . 
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Bethe-Salpctcr equation in which the kernel contains contributions, additional to 
those of the ladder truncation, whose form is constrained by the properties and 
consequences of the non-Abelian anomaly. The isotriplet solutions of this equation 
are precisely those of the ladder truncation: the non-Abelian anomaly does not 
alter the fact that the pion is QCD's Goldstone mode. However, no solution with 
the ideally-mixed flavour structures ud + du or ss exists. Instead, for the ground 
state pseudoscalar mesons the equation yields two solutions, at masses m no , m^, 
with m^ < m, )0 < m n > a , whose bound state amplitudes can be represented as: 
?/o - cose( m )F r°(k;P)-sm8( Vo )F s rfJh,Py, and r)' = sin 9 (ri' ) J 70 T° (k; P) + 
cos 0(r7 o )J rS r^, (k;P), where T° represents (uu + dd— 2ss)/y/& and JF S , (uu + dd + 
ss)/y/3. In this way one readily obtains a reasonable picture of the 77-77' complex, 
with m Vo « 3.3^1^, m n ' Q « 1.8m,, and #(770) = —10.5°, 0(r]' o ) = —31°. (NB. A 
mass-dependent mixing-angle is natural with the DSEsPSI Ideal mixing corresponds 
to 8 = —54.7°) In the chiral limit, rh u ,d = m s — 0, the model yields: 6(770) = 0° = 
0(t]' o ), and hence the 770 is purely U(3) octet while the r)' Q is purely U(3) singlet; and 
m-Ka = "tjjo = while m™ =0 = 0.9 m^. 

This model can be adapted to estimate the non-Abelian anomaly's effect on 
pseudoscalar meson radial excitations. Thus applied, it predicts m m « 1.08 m^, 
m^j w 1.22 m m and 6(771) = —35°, 6(r)[) = —43°, which matches expectation: the 
contribution to masses from QCD's nonzero topological susceptibility diminishes 
with increasing bound state mass. It follows that the 771's flavour content is 87% 
isosinglet, so that m Vl > m 7Tl . Moreover, the 77^ is 84% ss and hence the rainbow- 
ladder truncation is again a good approximation. This discussion suggests that 
77(1295) is the first radial excitation of 77(548) and, from line 3 of Tabled th & t 
7/(1475) might reasonably be interpreted as the first radial excitation of rf (958). 

The non-Abelian anomaly plays a nugatory role in the spectrum of mesons 
containing c-quarks. Hence the rainbow-ladder truncation should be realistic. T his 
may also be inferred from Ref. While the light masses were fixed elsewhere,^! 
the c-quark mass in Table was chosen to fit m^ng)- The values of the other 
quantities in Row 4 are a prediction. For model comparison (in GeV): Ref. |2U finds 
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Fig. 2. Evolution of meson masses with lj, where r a = X/u is the range of strong attraction in 
the model of Ref. 1171 . It is apparent that radial excitations are sensitive to the long-range part of 
the interaction, becoming less bound as the range of strong attraction decreases. 



f Vc(ls) = 0.213 and m Vc(2S) = 3.47, and Ref. 25 , /„ c(is) = 0.247. 

The model in Ref. |17j has a single parameter, u>. It may be used to define a 
length-scale r a = 1/lu: magnifying r a increases the range of strong attraction in the 
model. In this connection, Fig.[3]is significant. It illustrates the evolution of masses 
with uj and reveals plainly that the properties of the pion's first radial excitation 
are sensitive to the long-range part of the interaction; e.g., while a 30% increase in 
r a raises m wo by merely 3%, it reduces m 7ri by 14%. The sensitivity is explained by 
the evaluated charge radiusPI r- !Tl = 1.7 rv , which indicates that the pion's first 
radial excitation is a very extended object. This is materially important for many 
reasons. One instance is that a diameter d Wl = 2r Tri w 2.25 fm exceeds the spatial 
length, L < 2 fm, of lattices used in modern numerical simulations of lattice-QCD. 

Figure [5] also depicts the w-evolution of the mass of the kaon's first radial ex- 
citation and shows that mK„ =1 /m nri=1 w constant. We use this fact along with the 
7r(25)'s experimental mass to obtain mx(2S) : = m K, l=1 = 1-52 ± 0.12 GeV. This 
agrees with Ref. |24j . and suggests that the state identified as K(1A60) is the kaon's 
first radial excitation. 

The Dyson-Schwinger equations (DSEs) provide a unique perspective on non- 
perturbative aspects of QCD, e.g., they give: a straightforward understanding of 
the origin of constituent-quark masses; a veracious description of QCD's Goldstonc 
modes; and the essential connection between these phenomena. The existence of 
a nonperturbative systematic and symmetry-preserving truncation scheme enables 
the proof of exact results, which establish that chiral symmetry and the pattern 
of its breaking have a material impact even on the properties of hadrons with 
masses above 1 GeV. An appreciation of the essence of this truncation scheme en- 
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ables the formulation of a phenomenologically efficacious renormalisation-group- 
improved rainbow-ladder D SE model, corrections to which are quantifiable. The 
study of baryons is feasible!^ A discussion focusing on QCD's gauge sector may 
be followed from Ref. |27| . 
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